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Abstract 

We present the model-independent studies of non attractor inflation in the con¬ 
text of effective held theory (EFT) of inflation. Within the EFT approach two inde¬ 
pendent branches of non-attractor inflation solutions are discovered in which a near 
scale-invariant curvature perturbation power spectrum is generated from the interplay 
between the variation of sound speed and the second slow roll parameter r/. The first 
branch captures and extends the previously studied models of non-attractor inflation in 
which the curvature perturbation is not frozen on super-horizon scales and the single 
held non-Gaussianity consistency condition is violated. We present the general expres¬ 
sion for the amplitude of local-type non-Gaussianity in this branch. The second branch is 
new in which the curvature perturbation is frozen on super-horizon scales and the single 
held non-Gaussianity consistency condition does hold in the squeezed limit. Depending 
on the model parameters, the shape of bispectrum in this branch changes from an equi¬ 
lateral configuration to a folded configuration while the amplitude of non-Gaussianity 
is less than unity. 
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1 Introduction 


The basic predictions of models of inflation are well consistent with cosmological observa¬ 
tions DEI. The large scale curvature perturbations generated from quantum fluctuations 
during inflation are nearly scale-invariant, nearly adiabatic and nearly Gaussian as verified 
very accurately in recent cosmological observations. Despite the immense success of inflation 
as the leading paradigm for early universe and structure formation there is no unique theory 
of inflation. There are many models of inflation consistent with data. The simplest models 
of inflation are based on the dynamics of a scalar field rolling slowly over a near flat poten¬ 
tial. Even in the context of single field inflationary models, there are numerous scenarios of 
inflation. Therefore, it is an important question how far one can capture the most robust 
predictions of models of inflation without relying on particular realization of inflation. 

Effective Field Theory (EFT) of inflation (3] has provided a good answer to this question. 
In EFT language, all interactions compatible with the underlying symmetries are allowed. 
Then different models of inflation are realized depending on how one turn on particular 
interactions in the effective action governing the dynamics of the light field. So far the best 
studied set up of EFT of inflation is the single field models of inflation in which only one 
field is responsible for generating the curvature perturbations. However, the picture of EFT 
of inflation can be extended to more complicated setups such as multiple field scenarios [_4j. 

One interesting applications of EFT of inflation is its ability to classify different models of 
inflation based on their predictions for non-Gaussianity. Although there is no detection of non- 
Gaussianity from the Planck data [5], but non-Gaussianity continues to play important roles 
in constraining models of inflation at an age of precision cosmology mmmm- In particular, 
it is known that models of single field slow-roll inflation follow the Maldacena’s consistency 
condition mm in which the amplitude of non-Gaussianity /jvx in the squeezed limit is 
related to the spectral index of curvature perturbation power spectrum n s via f^ L ~n a -10 
This has led to the conclusion that a detection of local type non-Gaussianity at the order 
Inl = 0(1) w iH rule out all single field models of inflation. However, it was shown in [T5] 
that this general statement does not hold in models of non-attractor inflation in which the 
slow-roll parameter e = — -J^-, measuring the variation of the Hubble parameter H, falls off 
exponentially during inflation. As a result, the would-be decaying mode of the curvature 
perturbations actually grows during the non-attractor phase and the curvature perturbations 
1Z is not frozen on super-horizon scales. Therefore, the single field non-Gaussianity consistency 
condition is violated and f nl > 0( 1) can be generated during the non-attractor phase jT3l 

mm- 

So far there are few known examples of non-attractor scenarios. The simplest model is 
presented in [13], see also mm, in which the scalar field is rolling in a constant potential 

1 It is worth mentioning that, as argued in |12j . the freely falling observers in Fermi normal coordinates 
can not observe this small amount of non-Gaussianity as bispectrum in these sets of coordinates is cancelled 
up to corrections of 0(k\/k |) in which and ks are the long mode and the short mode defined in squeezed 
configuration. 
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with sound speed c s — 1. A different model of non-attractor inflation was presented in 
m 153 i n the context of P(X, 0) setup of K-inflation [TS] IT? ]. see also [2D], with a constant 
sound speed c s . With the motivation of EFT of inflation, it is a natural question to ask 
what are the generic predictions of non-attractor models of inflation without relying on a 
particular realization of non-attractor scenario? The goal of this paper is to address this 
question. As we shall see, the non-attractor scenarios can easily be embedded within the 
general context of EFT of inflation. Equipped with the power of EFT of inflation, we present 
the general predictions of the non-attractor scenarios for the spectral index and the shapes of 
non-Gaussianities. We demonstrate that there are two independent branches of non-attractor 
models which can generate a nearly scale-invariant power spectrum. In one branch we capture 
and extend the previous results obtained in [T3j [14i, |T5] . Very interestingly, the other branch 
of non-attractor scenarios is new which escaped the previous model-dependent studies. 

2 The non-attractor inflation in EFT setup 

In this section we present our setup of non-attractor inflation in the context of EFT of 
inflation. Let us first review the EFT of inflation very briefly. An extensive discussion can be 
found in [3]. In what follows, we use the conventions and the methods employed in [21] . 

In a quasi-de-Sitter background with a time-dependent scalar held 0(f), the four-dimensional 
diffeomorphsim invariance is spontaneously broken to a three-dimensional spatial diffeo- 
morphsim invariance determined by 0(f) = constant. In the so-called unitary (comoving) 
gauge, one is allowed to write down all terms in the action which is consistent with the 
three-dimensional diffeomorphsim invariance while yielding the known backgroimd for the 
Hubble expansion rate H(t) and its derivative H(t). Alternatively, one can change this pic¬ 
ture by restoring the explicit full four-dimensional diffeomorphsim invariance by introducing 
a scalar held fluctuations, the Goldstone boson associated with the breaking of the 

time diffeomorphsim invariance. 

In the presence of the perturbations, the metric is given by 

ds 2 = —N 2 dt 2 + gij(dx l + N l dt)(dx j + N j dt) (1) 

in which V(x M ) and N l (x lL ) respectively represent the lapse and the shift functions and gij 
represents the spatial metric. 

The full action, after restoring the four-dimensional diffeomorphsim invariance by the 
introduction of the Goldstone boson ^(a^), is Slotai = Setter + Seh in which the matter part 
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of the action is 


5, 


matter 


= Jd A Xy/^g —MpH(t + 7 r) (l + tt - A^<9;7t) 2 - g 13 

-Ml (^3 H\t + vr) + H{t + tt)) + 

M(t 2~ ( l + * ~ N ' d ~ 9 ij d l 7rd J 7i - l) + 

c 3 (t + tt) M(t + 7r) 4 / 1 


( 2 ) 
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^■^2 ( X + * ~ N'diir ) 2 - g^diTrdjTT - 1^ + .. 


In this picture the first two lines of Eq. (J2D are fixed by the given form of the FKW background 
while the term containing M(t + n) in the third line controls the sound speed of scalar 
perturbations, c s |3]. More specifically, the relation between M[t ) and c s is given by 



Note that both operators in the third and fourth lines of Eq. (J2J) containing M(t +7r) originates 
from the term g 00 +1 in the unitary gauge (comoving gauge) and c 3 {t + 7r) is an undetermined 
coupling. As understood from the logics of EFT, there are many terms which are not written 
in the above action as they are not important in our limit of interest, the decoupling limit, in 
which the gravitational back-reactions are discarded and when we are well below the strong 
coupling limit. 

In addition, the gravitational part of the action from the Einstein Hilbert term is 

S EH = \m 2 p J d A x yf^g R = X -Ml jd 3 x dt ^ [NR® + ^(H 3 E l} - E\ 2 )] , (4) 

in which Mp is the reduced Planck mass, E,j is related to the extrinsic curvature K l3 associated 
with the three-dimensional hypersurface 0(f) = constant via E i3 = NK i3 and R® represents 
the three-dimensional Ricci scalar associated with the metric g l3 . 

As discussed in |3j the great advantage of EFT of inflation is in the decoupling limit in 
which one can safely neglect the gravitational back-reactions so the observational predictions 
such as the tilt of scalar perturbation or the non-Gaussianity parameter f NL are determined 
to leading order by the matter fluctuations. Technically, this corresponds to setting N = 1 
and N l = 0 in the action and neglecting the perturbations from the gravitational action Eq. 
PD compared to perturbations from the matter action. Having this said, in the appendix 
we integrate out the non-dynamical variables N and N l and demonstrate the validity of the 
decoupling limit for our analysis. 

After presenting the setup of EFT, here we present our definition of non-attractor scenar¬ 
ios. For this purpose, it is convenient to define the “slow-roll” parameters associated with the 
variations of H(t) via 


H 




H 2 



(5) 
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Our general definition of the non-attractor scenarios is that during non-attractor phase the 
slow-roll parameter e falls off exponentially such that the background becomes more and more 
like a dS background as inflation proceeds. Consequently, the second slow-roll parameter 77 
is not small. Indeed, we expect 77 ~ 0(—1) during the non-attractor phase. Not that the 
requirement that e falls off implies that 77 < 0. The exact value of 77 is determined from n s . For 
example, in non-attractor models studied in [13] [IT, 15] in which the sound speed is constant, 
n s is given by n s — 1 = 6+77 so requiring a near scale-invariant power spectrum to be consistent 
with observations implies 77 ~ — 6 . In these models, with 77 ~ — 6 , we find that e falls off like 
l/a(t) 6 . As we shall see, during this non-attractor phase the curvature perturbation 1Z is 
not frozen on super-horizon scales and it grows like a(f) 3 . This is the main reason why the 
celebrated Maldacena’s consistency condition is violated in non-attractor model. We stress 
that the non-attractor phase can not extend for a long period, as e becomes very small while 
1Z grows exponentially and the system becomes non-perturbative. In order to prevent this to 
happen, the non-attractor phase has to be followed by an attractor phase in which e becomes 
nearly time-independent and 1Z saturates on super-horizon scales. 

Motivated with the generality of EFT, here we extend the non-attractor models to case in 
which the sound speed and the operator c^it) can also vary during the non-attractor phase. 
For this purpose, let us define the slow-roll parameters in the matter sector via 


Hc s 



( 6 ) 


In principle the parameters s or d 3 can take arbitrary values as long the power spectrum 
becomes nearly scale-invariant. 

With these discussions one may ask what conditions should be imposed on the variations of 
the secondary slow-roll parameters 77 , s and £ 3 . We take the simple assumption that these sec¬ 
ondary slow-roll parameters are nearly constant and ?) ~ s ~ £3 ~ 0(e). Of course, this may 
not be the case in general and these secondary slow-roll parameters may have rapid variations 
like e, c s or c 3 , but as we shall see, even the simple assumptions of allowing ( 77 !, |s|, |< 5 3 | 0 ( 1 ) 

while neglecting their evolution describes large enough class of non-attractor models. 

This discussion summarizes our definition of non-attractor phase without relying on any 
particular model. Our goal is to understand the properties of non-attractor models in a 
model-independent way using the effective held theory (EFT) method of inflation |3|. As we 
shall see our model-independent analysis reproduces the results obtained in [T3( [IT], 15] in 
particular limits. 

Now the question arises how the non-attractor models can be embedded in EFT approach 
of inflation. Interestingly, this is vary simple. During the non-attractor phase e falls off 
rapidly so the the main places where the effects of non-attractor phase enter are the first two 
lines of action ([2]) containing H(t + n) and H 2 (t + 7 r). Therefore, we need to expand H(t-\- n) 
and H 2 (t + 7 r) to first order of e. Note that the presence of 77 ~ 0(— 1) will play a crucial role 
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here. We have 


and 


H = -eH 2 -2eHH 

= -r]eH 3 + 2e 2 H 3 =-r]eH 3 + 0(e 2 ), 


— = — fjaH 3 - r] 2 eH 4 + 7rje 2 H 4 - 6 e 3 H 4 = -ifeH 4 + 0(e 2 ) 
at 


(7) 


( 8 ) 


Note that in the last relation we have neglected the evolution of r] because = 0(e) as 
assumed above. 

Also, 

d 2 


dt 2 


H = -rfeH 3 + 0(e 2 ). 


(9) 


Combining these expansions for ^ and to leading order for H 2 (t + 7 r) we obtain 


H 2 (t + vr) = H 2 (t) + 2HHtt + - (h 2 + 2HH] i r 2 


1 ( 


6 


HH + 2 H 


dH 

dt 


7T 3 + . . . 


= H 2 (t) - 2eH 3 7T - r)eH 4 n 2 - -r] 2 eH 5 ir 3 + ... . 

3 

Similarly, for H(t + n) to leading order we obtain 

H(t + 7r) = —eH 2 - r]eH 3 iT - -erj 2 H\ 2 - Vetf V + ... . 

2 6 


( 10 ) 


( 11 ) 


We stress that in obtaining the above expansions, we have kept terms leading in e because e 
falls off rapidly during the non-attractor phase. In addition we have neglected the variations 
of 77 as we assumed that p 0(e). 

In addition to contributions from H 2 (t + ir) and H(t + n) from the background dynamics, 
we also have the contributions from the rapid variations of the sound speed c s (t + 7 r) and the 
coupling c 3 (t + 7r). These contributions should also be included in quadratic and cubic actions. 
As we discussed before, the operator M(t) controls the sound speed of scalar fluctuations as 
given in Eq. ([3]). As a result 


1 - Hu + 77 I 1 -5- ) Htt 


M 4 (t + tt) = —~eH 2 Mp 


Similarly, for the operator c 3 (t) we have 

c 3 (t + 7 r) = c 3 (t) (1 + 5 3 Htt) + ... 


+ ... 


( 12 ) 


(13) 


Equipped with the expansion of the background parameters H(t + n) and H(t + n) as given 
in Eqs. (flOl) and (ITT]) and the matter sector parameters Eqs. (fl 2 |) and (fl3|) we can obtain 
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the quadratic and cubic action for n in the matter sector along with the action from the 
gravitational sector from Eq. (HI) . 

As we mentioned before, the advantage of EFT of inflation is in decoupling limit in which 
one can neglects the gravitational back-reactions. This is even more justified in the non¬ 
attractor regime in which e falls off rapidly and the errors in neglecting the gravitational 
back-reactions which are controlled by higher powers of e are quite negligible. As we see from 
the detail analysis in the Appendix, we can safely go to decoupling limit in which we can 
ignore the contribution from Eq. (3j) and set N — l,N l — 0 in the matter sector. In this 
limit, the quadratic and the cubic actions generated from Eq. ([2]) respectively are 



The above actions are written in terms of n. However, we are interested in power spectrum 
and bispectrum of comoving curvature perturbation 1Z. The relation between n and 77, to 
quadratic order which is necessary for the bispectrum analysis, have been worked out in 
' IT) . 21 ] yielding 

TZ = —Htt + H'k'K + —it 1 + 0(7T 3 ). (16) 

We have specifically checked that the above relation does hold for the non-attractor setups 
too. 

Finally, to calculate the quadratic and cubic actions and the follow up in-in analysis we 
need to know the functional form of the quantities e(t),c s (t ) and c^{t). Using the definition 
of rj, s and 83 given in Eqs. (J5]) and ([ 6 ]), and employing the assumption that 77 , s and £3 do 
not vary during the non-attractor phase we have 



in which r is the conformal time related to cosmic time via dr = dt/a(t ) and r e represents 
the time of end of non-attractor phase. Finally, the scale factor in terms of conformal time is 

given by a = a e . 

As we shall see, a combination of 77 and s determine the tilt of curvature perturbation 
power spectrum. Therefore, neither rj nor s is fixed individually. Having this said, we can 
derive an upper bound on the value of s as follows. The universal definition of inflation is 
that a > 0 so during this phase the comoving Hubble radius 1 /aH falls off rapidly. However, 
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in the inflationary model with a non-trivial sound speed, we encounter the “sound horizon” 
associated with the scalar perturbations. Correspondingly, the comoving sound horizon is 
given by c s /aH. In order for the perturbation with the comoving wave number k to be 
sub-horizon during early stage of inflation and then to leave the comoving sound horizon 
during the subsequent stage of inflation, we require that the comoving sound horizon to be 
a decaying function during inflation. Using Eq. (113, the condition for the comoving sound 
horizon c s /aH to be a decaying function as inflation proceeds is s < 1 which will be imposed 
in our analysis below. 


3 Power spectrum 


In this section we study the predictions for curvature perturbations power spectrum. To 
leading order 1Z = —Hu and the quadratic action for curvature perturbations power spectrum 

is 


So — - I cPxdrz 2 


n 


/2 




(18) 


in which a prime indicates the derivative with respect to conformal time, d,7Z represents the 
spatial derivative of 7 Z and the parameter z is defined via 

2 ea 2 


z 2 = 


.2 - M l' 


(19) 


Interestingly the action (fTHjl coincides with action obtained in particular P(X , 0) model of 
non-attractor inflation studied in mm- In the decoupling limit in which we neglect the 
gravitational action (J4j) and, as long as we neglect the variation of rj, the second order action 
for all non-attractor scenarios is uniquely given by Eq. (TT8|) . The only relevant parameter is 
the sound speed of perturbations c s (t), which is determined by the operator M 4 (f) given in 
Eq. (H2D- 

To quantize the system and to calculate the power spectrum, it is conventional to define 
the canonically normalized field v related to 1Z via v = zTZ. After some integrating by parts, 
the action for the canonically normalized field v is given by 


So = - 


/ d 3 xdr 

iy') 2 - c 2 Mv) 2 + -v 2 

1 

z 


( 20 ) 


in which the quantity —represents the effective time-dependent mass of the canonically 
normalized field. More specifically, to leading order of e we have 


z" 

z 


r z 


2 + ~(r] - 2s)(rj - 2s + 6) 


( 21 ) 


The equation of motion for the canonically normalized field in the Fourier space with the 
momentum number k is 


v'l + I c 2 k 2 - 


V 2 - 


T z 


1 

4 


Vk = 0 , 


( 22 ) 








in which a prime indicates the derivative with respect to the conformal time and the index v 
is defined via 

_ r\ — 2s + 3 


We present the exact solution of Vk(r) in Eq. (1221) for the form of c s given in Eq. (fTT|) . 
However, before presenting the general solution for Vk(r) let us pause and check whether or 
not the notion of propagating wave and the initial adiabatic mode is well-defined in our system 
when |s| is not small and c s (r) evolves rapidly. In order to have a well-defined propagating 
mode deep inside the horizon, we impose the WKB approximation in which \d s k/{c s k) 2 | <C 1. 
For this to apply we require 



c s k\r\ |s 


(24) 


In this limit, one can neglect the second term in the big bracket in Eq. (12211 and a WKB 
solution for modes propagating deep inside the horizon is applicable. The corresponding 
WKB solution is 

v(k,r ) ~ 


Note that the only change compared to standard case in which one can neglect the evolution 
of c s is the appearance of the additional factor 1 — s in the denominator of the argument. As 
we discussed before, in order to have a decreasing function of comoving sound horizon, we 
require 1 — s > 0. 

Happily Eq. (122|) can be solved analytically. The solution which matches the WKB 
solution for modes deep inside the horizon Eq. (EZ5D is 

<*> 

in which Hj}\x) is the Hankel function of type one, x = —c s (r e )/cr and 


y /2 kc s 

1 

y/ 2 kc. 


exp 


: exp 


—ik J c s (k,T)dr 
■ik c s (t)t 


(■ c s k\r\ » \s\) . 


(25) 


is 


/i = 


C', = 


7r 


D i(M+l/2) 7 r/2 


4fcc*(l 


(27) 


We are interested in the power spectrum of the curvature perturbation TZ = v on 

super-horizon scales in which x —> 0. Using the following formula for the small argument of 
Hankel function 


H£\x) = e -<0*-MW2 H$(x) ~ 


-xm f2\ 

IT \X J 


ImI 


(x —> 0) 


(28) 
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the curvature perturbation on super-horizon scales is obtained to be 


K(t ->■ 0 ) 


2H-§r(|/i|)(l-s)W-3 H (t 


V 


7 re P c« 


U3 

o IV 




T P 




- c s kr) 


l-H 


(29) 


in which e e and c se represents the value of the corresponding quantities at the end of non¬ 
attractor phase r = r e and T is the Gamma function. 

Correspondingly, the power spectrum at the end of non-attractor is obtained to be 


Vn(r e ) 


k 3 

27T 2 


K(r e ) 



(30) 


in which k e = a{r e )H / c se represents the mode which leaves the sound horizon at the end of 
non-attractor phase and the amplitude A (COBE normalization) is 


= r(|^|) 2 (! - 1 /_H_y 

5 , 324 - 21^1 (eCse \Mp) 


(31) 


The spectral index n s from the power spectrum (130|) is obtained to be n s 
Depending on the sign of /i we have 

- 1 = 3-2|/i|. 

n s 1 - (h > 0 ) 

(32) 

1 — s 

or 


n,-i = 6 +’'- 5s m < 0). 

1 — s 

(33) 


The branch represented by Eq. (1321) reproduces the results in the known slow-roll limit, the 
K-inflation model [T9J, [35, !36], in which |s|, |r;| -C 1 and n s — 1 ~ —rj — s — 2e . Note that 
we have already discarded e in our analysis assuming that e —* 0 so the 0(e) discrepancy 
with the K-inflation results is expected. Furthermore, the branch represented by Eq. (Holl is 
a generalization of the previously studied models of non-attractor with s = 0 H3HHII5] in 
which 77 ~ —6 to get a nearly scale-invariant power spectrum. Note that the branch r/ ~ — s 
represented by Eq. (l32l) is new which was not noticed in the previous model-based analysis 
of PJ, I3| dSJ . This demonstrates the power of EFT as a platform to study inflation model 
independently. We comment that the condition 77 = — s to obtain a scale-invariant power 
spectrum may be interpreted from the fact that Vn oc 1 /ec s . So, one may expect that a 
variation of c s (t ) should be balanced by the variation of e(t) such that the power spectrum 
remains scale invariant when 77 = — s. 

Equipped with the powerful EFT method we see that there are different options to get 
a nearly scale-invariant power spectrum with a mild red-tilted power, n s ~ 0.96, depending 
on whether 77 > 0 or /i < 0. Note that our starting assumption was that e falls off rapidly 
so we choose the sign of 77 to be negative. As for s, as we discussed before, we only require 
s < 1 in order to get a decreasing function of comoving sound horizon during inflation so 
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the sign of s is undetermined. A positive (negative) s indicates a growing (decaying) sound 
speed during inflation. In the case s > 0 the sound speed can cross unity and we enter the 
superluminal regime. In order for this not to happen, one concludes that the duration of the 
non-attractor phase to be finite, say few e-folds such that |siV e | < 1 in which N e represents 
the duration of non-attractor phase. On the other hand, if s < 0 then the sound speed falls off 
as inflation proceed. If c s becomes arbitrarily small, then one encounters the strong coupling 
regime limit with very large non-Gaussianity. Therefore, in order to avoid the strong coupling 
limit, we demand that the non-attractor phase terminates before c s becomes very small, say 
c s > 0.003 (3j. However, there is stronger upper bound on c s from the Planck’s constrains on 
equilateral-type non-Gaussianity, requiring c s > 0.02 |[5j. 

Having obtained the two possible branches to get a nearly scale-invariant power spectrum 
as give by Eqs. (I52|) and (j33|) . it is instructive to look at the evolution of 1Z on super-horizon 
scales. From Eq. (123|) . and noting that \fi\ ~ | to get a scale-invariant power spectrum, we 
obtain 

7 Z= - H1Z (kr —>■ 0). (34) 

Interestingly, for the branch r) = —s, the curvature perturbation is frozen on super-horizon 
scales. This is similar to the conventional models of single held slow-roll inflation. In addition, 
as discussed in mm, the constancy of 1Z is the essential reason for the validity of the single 
held non-Gaussianity consistency condition. Therefore, in this branch we expect that the 
single held non-Gaussianity consistency condition to hold and = O(e). We verify this 
conclusion explicitly in next Section. On the other hand, for the branch 77 = —6 + 5s, we 
obtain 1Z = 3(1 — s)H1Z on super-horizon scales. This is in agreement with the results 
obtained in the particular P(X , 0) model studied in [T3l[Til[T5] in which s = 0 and 1Z = 3HTZ 
on super-horizon scales. As a result, since 1Z is not frozen on super-horizon scales, we expect 
that the non-Gaussianity consistency condition to be violated in this branch as we shall see 
explicitly from our in-in analysis. Note that since s < 1 we conclude that 1Z > 0 and the 
curvature perturbations grows exponentially during the non-attractor phase in this branch. 
This indicates an instability so we have to terminate the non-attractor phase followed by an 
attractor phase P3J- 

As mentioned above, the branch 77 = —6 +5s is an extension of the previous non-attractor 
models based on P(X,(p ) action as studied in [13], (14, IS]. Therefore, it is an interesting 
exercise to realize a specific model for the branch 77 = —s such as in P(X, 0) setup. In 
principle one could have discovered the branch 77 = —s from examples built in P(X, 0) setup. 
As experienced in H3 this requires careful inverse engineering. Having this said, we find it 
interesting that the new branch 77 = — s emerges naturally within the general context of EFT 
of inflation. 
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4 Bispectrum 


Having studied the power spectrum and the criteria to obtain a nearly scale-invariant power 
spectrum, here we study the predictions of the general non-attractor models for non-Gaussianity. 
Here the discussions become more interesting. First, EFT provides a framework in which dif¬ 
ferent inflationary models can be classified based on their predictions for the amplitude and 
the shape of non-Gaussianity. Second, the particular non-attractor models presented origi¬ 
nally in [131 EES IS] are the very few known examples of single field inflation which violate the 
non-Gaussianity consistency condition, for more discussions on non-Gaussianity consistency 

condition see mmmmmmmmmmmm- Therefore, it is natural to ask if this 

violation of non-Gaussianity consistency condition is a generic feature of the non-attractor 
inflationary systems. 

To calculate the bispectrum of TZ we have to take into account the non-linear relations 
between TZ and tt given in Eq. ([T 6 lh More specifically, the bispectrum of 7 Z is related to the 
bispectrum of n via 


^(k 1 )7^(k 2 )^(k 2 )^ = -H' 3 ^ 7 r(k 1 ) 7 r(k 2 ) 7 r(k 3 )^) + H' 3 ^ 7 r(k 1 ) 7 r(k 2 ) ( 7 T 7 t) (k 3 )^ + 2c.p. (35) 
The last term above involves convolution integrals, yielding 


ir(ki)7r(k 2 ) (mr) (k 3 )) = — [|R(k,)| 2 R(k 2 )R(k 2 )‘ + Iftlfe^RfkORfk, 


(36) 


in which ()' indicates that we have absorbed the common factor (27t) 3 <5 ( - 3 ' > (k 3 + k 2 + k 3 ). Now, 
using Eq. (1341) to express TZ in terms of TZ, the contribution in bispectrum from the non-linear 
relation between n and TZ is 


H 3 


^ 7 Tki 7 Tk 2 3 ^ + 2c.p. — — (rj + s) P 7 e(ki)Pft(k 2 ) + P^k^P^^) + P-£.(k 2 )P£.(k 3 ) (37) 


Plugging Eq. (1371) in Eq. (|35lh the relation between the bispectrum of tt and TZ is 

Jp(k 1 )P(k 2 )P(k 2 )) / = ^7r(k 1 )7r(k 2 )7r(k 3 ) 

- (77 + s) \p n (k 1 )P'i Z (k 2 ) + P7e(ki)P7e(k 3 ) + P n (k 2 )PnO^s) (38) 


in which it is understood that all quantities are calculated at the time of end of non-attractor 
phase r e . The power spectrum P^(k) at r e is 


P 1 z (k,re) = \TZ(r e )\ 2 


(1 - s) 2 H 2 
4e e c e k 3 Mp 


(39) 


To calculate the first term above, ( 7 r(k 1 ) 7 r(k 2 ) 7 r(k 3 )y, we need the cubic action which is 
given in Eq. fjT5]l . Specifically, using the standard in-in formalism [371 [ 8 ], [38] the bispectrum 
at the end of non-attractor phase is given by 


7r(ki,T e )7r(k 2 ,r e )7r(k 2 ,T e )^ = i / dr^7r A . 1 (r e )7r fc2 (r e )7r fc3 (r e )£i(r)^ + c.c. 


(40) 
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in which c.c. stands for complex conjugation and A, represents either of the four Lagrangian 
terms in Eq. (fT5|) : namely 


Ci 

c 2 

c 3 

u 


H : 


Xv - 2s )« 2 


eMpTin 


/ 2 


-r]H' i a 2 Mpe'K(diTTY, 


-aH 2 


2 

i + jCs 


1 - \ ) eVV 3 , 



(41) 

(42) 

(43) 

(44) 


Since the details of the analysis and the shape of non-Gaussianity are distinctly different 
for the two allowed branches of scenarios r) = —6 + 5s and r] = —s, we study each branch 
separately. 


4.1 The branch rj = —s 


In this case the wave function of 1Z is given by 


TZ(k, t) 


2 aM p 


7 TT 


2 e(l 




(45) 


in which we have neglected the corrections at the order 0(1 — n s ) in wave function and set 
H = |. This is justified since we are interested in generating f NL ~ 0(1) and can discard the 
sub-leading corrections to f nl at the order n s — 1 or e. 

Happily the in-in integrals can be performed analytically. We present the result term by 
term. For the contribution from C\ we have 


^(k!)^(k 2 )^(k 3 


3 H 4 


s s 


l) 4 


X 


\m%K^k\k\k\clel 

cos(27ts) (^kik 2 kY) (k\k 2 + k\k 3 + k 2 k 3 ) + K{k\k1 + k\kl + k^k^) 


(46) 


in which we have defined K = k 3 + k 2 + k 3 . 
The contribution from C 2 yields 


^(k 1 )7e(k 2 )7e(k 3 ))^ = 

x 


H\k\ + k\ + k\) s(s-l ) 4 
32M|,A' 2 elclk\k\k\ 

cos(27t s)k±k 2 k 3 — K(k\ + k 2 + k 2 + k\k 2 + k\k 3 + k 2 k 3 ) 


Finally, C 3 and £4 make no contributions (TZ 2 )^ = (A 3 )/^ = 0. 

We also have to include the contributions from the non-linear relation between n and 
1Z. Interestingly, in the current case in which 77 + s = 0, there is no contribution from the 
non-linear relation between n and 1Z as can be seen from Eq. (137]) . As a result, the total 
contribution to bispectrum is from C\ and £ 2 . 
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To simplify the analysis, it is useful to define 


K = ki + k 2 + k 3 , P 2 = fcifc 2 + hk 3 + k 2 k 3 , Q 3 = kik 2 k 3 . (48) 


With these definition, the total bispectrum is 


(KXm = 


S [S 


1) 4 P 4 


32Mp6 2 c 2 K 2 Q 9 L 


Q 3 (/W + 4P 2 ) cos(2vrs) + K (4P 4 - K 4 + 3P 2 P 2 - 12 Q 3 K ) (49) 


The above bispectrum has a non-trivial shape. However, from Eq. (134|) . we know that with 
r/ + s = 0, the curvature perturbation 1Z is frozen on super-horizon scales. Therefore, the non- 
Gaussianity consistency condition should be satisfied in the squeezed limit mm in which 
/wl 0(1 (1 ~ n s) 0- Therefore, it is constructive to look at the shape of bispectrum Eq. (149|) 
in the squeezed limit k 3 k\ ~ k 2 - In this limit K ~ 2k\, P 2 ~ /c 2 and Q ~ k 3 k\ — > 0 and 
Eq. (1491) yields 


(K 3 )'Z* ^ s^P K (h)P K (k 3 ) ^ 0. 


(50) 


so the bispectrum vanishes in the squeezed limit, in agreement with the Maldacena’s consis¬ 
tency condition. Note that we have neglected the gravitational back-reactions and also set 
n s = 1 in the bispectrum analysis so we can not recover 0(e) corrections in bispectrum in the 
squeezed limit. 

It is also instructive to look at the magnitude of non-Gaussianity in the equilateral limit 
k\ — k 2 — k 3 . In this limit we obtain 

(P 3 )' eq = ^-s cos(27rs)P^(fci). (51) 

6 

A convenient way to parameterize the amplitude of non-Gaussianity is via the parameter f^L 
defined via 


6 (K(k 1 )K(k 2 )K(k 3 ))' 

5 fm ~ P^iMWM) + 2c.p. ' 


(52) 


Using this definition, the non-Gaussianity parameter in the equilateral configuration f^ L is 
obtained to be 

35 

Inl = Y^ scos ( 27rs )- ( 53 ) 

Knowing that s < 1 in order to have a decaying sound horizon during inflation, the above 
equation indicates that the amplitude of f^ L is less than unity so it is consistent with obser¬ 
vational constraints. 

Now let us look at the shape of bispectrum in general case. Following the convention of 
Chen [8], the dimensionless shape function S(ki, k 2 , k 3 ) is defined by 
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Figure 1: The amplitude of the shape functions, (S'!, given in Eq. (l54|h Left panel is for 
s = 0.9 so the term containing cos(27ts) dominates yielding a nearly equilateral shape. Right 
panel is for s = 0.3 and the shape is close to the folded shape. Our convention is that Xi = 
for i = 2, 3. 


yielding 


S(k 1 ,k 2 ,k 3 ) 



4P 4 K 3 | 3KP 2 

kq>~q j + ~q t ~ 



(54) 


Note the curious form of the numerical factors. As can be seen from our shape plots, the 
numerical factors in the second big bracket above are such that the bispectrum does not 
peak in the local shape. We note that the term containing cos(27ts) generates the equilateral 
shape while the terms in the second bracket generate the folded shape. Therefore, depending 
on the value of the parameter s, the shape of non-Gaussianity changes from the equilateral 
shape to folded shape. The equilateral shape dominates when | cos(27ts) | is large, i.e. when 
s ~ 0 or s ~ 1, while the folded shape dominates when cos(27ts) becomes small when s ~ 
This is similar to the phenomena observed in model of quasi single field inflation [32| ;3T] , see 
also [33j [34], in which the shape function evolves from the equilateral shape to local shape 
depending on the mass of the semi-heavy iso-curvaton field. Finally, the amplitude of non- 
Gaussianity is less than unity for both equilateral and folded shapes so the predictions of the 
model for non-Gaussianity are well-consistent with observational bounds. We stress again 
that the shape does not peak for local-type configuration which is supported from the fact 
that the non-Gaussianity consistency condition is satisfied in the squeezed limit as discussed 
in Eq. (15Ujh 
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4.2 The branch rj = —6 + 5s 

Here we calculate the bispectrum for the case rj = — 6 + 5s. Only C\ and £ 3 yield non-zero 
contributions for in-in integrals and we have 


„3U r,-2s( s -l)VCfr; + *j + frf 
' lc ' 32 


(55) 


and 


Wc, = 


X 


(s - 1 ) 6 H 4 k\ + kl + kl 
e l c l e Mp k\klkl 

■ 9 6(1 — Cg) + 6sc 2 e — 4s 9c 3e (6 — d 3 )(l — c 2 se ) + 6sc^ e — 4s~ 
.32 (s — 1) (4s — 6) 8 (d 3 + 6s — 6) (4s + <f 3 — 6) . 


(56) 


in which c 3e represents the value of c 3 at the time r = r e . 

The contribution from the non-linear relation between it and 71, (TZ 3 ) n ^Ti, from Eq. (137|) 
is 






b,o Ld b,o 

^ 1 ^ 2^3 


(57) 


Adding the contributions from (JZ i )c 1 -, (7l 3 )c 3 and (7the final bispectrum is obtained 
to be 


(K 3 Y to tai = P(ki)P{k 2 ) + P(h)P(k 3 ) + P(k 2 )P(k 3 ) 


T 


(58) 


in which the function T is defined via 

r _ ~( r l + 4s) , 2 r 27(l-cg e )+9s(3cg e -2) (6 - d 3 )(l - c 2 se ) + s(6cg e - 4) ~ 

2 L (s — 1) (4s — 6) 3e (<5 3 + 6s — 6) (4s + S 3 — 6) - 


From Eq. (1581) we see that the bispectrum has the exact local shape and no other shapes are 
generated. 

Now let us calculate the non-Gaussianity parameter /jvl in the squeezed limit h < h ~ 
k 2 . Using the general definition given in Eq. ([52]) . f^L in the squeezed limit is obtained to 
be ffNL = 2jr yielding 

jf%L = -fa + 4s) (59) 

54(1 — c 2 e ) + 18s(3c^ e — 2) (6 — d 3 )(l — c 2 e ) + s(6c^ e — 4) " 

(s — 1) (4s — 6) 3e (d 3 + 6s — 6) (4s + d 3 — 6) - 

As an example consider the P(X, 0) model presented in [13] and [15] with s = S 3 = 0 and 
77 = — 6 . In this case we obtain 


+ ( s — l) 2 


J 7 = 6 + 3(1 — c^)(2c 3 + 3). 


(60) 
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Furthermore, if we assume C 3 = — | + ^2 which is valid in P(X, (ft) model studied in [14] and 
[15 ] then T = and 


Inl = 


5(1 +eg) 
4c? 


(61) 


in exact agreement with the results of ( 14] and [15]. Finally, if we further assume c s = 1, as 
in model studied in [13], then f^L — § which is well-consistent with the Planck constraints 
on the amplitude of local-type non-Gaussianity [5]. 


5 Summary and Discussions 

In this paper we have presented a model-independent study of non-attractor inflation in 
the context of EFT of inflation. The goal was to find the generic predictions of the non¬ 
attractor scenarios for the power spectrum and bispectrum. We believe that it is important to 
understand the general predictions of the non-attractor scenarios in a model-independent way. 
This is mainly because the non-attractor scenarios provide the very few examples in which 
the Maldacena’s non-Gaussianity consistency condition for single held inflation is violated. 
As a result, a detection of local-type non-Gaussianity with f NL ~ 0(1) does not necessarily 
rule out all single held models of inflation. More precisely, it rules out all single held models 
of inflation which has reached the attractor phase. 

We have studied a large class of non-attractor model in the context of EFT of inflation. 
Our only assumptions were that 77 and s do not evolve significantly during the non-attractor 
phase. These assumptions were imposed mainly to keep the analysis under analytical control. 
In principle one can go beyond these simplifications and allow the situations in which both 
77 and s have significant evolution during the non-attractor phase. Having this said, the 
setup studied here in which \r]\, |s| ~ 0 ( 1 ) while 77 , s ~ 0 represents large enough class of 
non-attractor scenarios. 

We have obtained two branches of non-attractor scenarios which can generate a near 
scale-invariant power spectrum determined by 77 = —6 + 5s or 77 = — s. The first branch 
encompasses the previously known models of non-attractor scenarios in which s = 0 and 
77 = —6 [T51 [T4l [ 15 ] . However, the branch 77 = —s is new which was not noticed previously. 
This demonstrates the power of EFT of inflation as a very helpful setup to study the paradigm 
of inflation model-independently. 

We have studied the predictions of the above branches of non-attractor scenarios for the 
bispectrum. We have shown that the bispectrum in the new branch 77 = —s has a non¬ 
trivial shape. Depending on the value of the parameter s, the bispectrum can change from 
the equilateral shape dominated when |cos(27ts)| —>■ 1, to folded shape dominated when 
cos(27ts) —* 0. For both configurations, the amplitude of non-Gaussianity is less than unity 
which is consistent with the observational constraints. In addition, in this branch 1Z is frozen 
on super-horizon scales and as a result the non-Gaussianity consistency condition is expected 
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to hold in this branch. We have checked explicitly that this is indeed the case. On the 
other hand, for the branch r] = —6 + 5s, only the local-type non-Gaussianity is generated. 
The curvature perturbation is not frozen on super-horizon scales and the non-Gaussianity 
consistency condition is violated as previously noticed in particular examples studied in [ T3l 
□3115]. Depending on the value of model parameters such as s, c se , 5 and c 3e a large local- 
type non-Gaussianity is generated. Combined with the Planck’s constraints on local-type 
non-Gaussianity this can be used to constrain the model parameters. 

An interesting conclusion of this model-independent study was that not all single field 
non-attractor models can violate the non-Gaussianity consistency condition. Only in non¬ 
attractor models in which the curvature perturbation is not frozen on super-horizon scales 
the single field non-Gaussianity consistency condition is violated. 

Acknowledgments: We would like to thank A. A. Abolhasani, X. Chen, R. Emarni, E. 
Komatsu and M. H. Namjoo for useful discussions and correspondences. 

A The quadratic and cubic actions 

Here we present some details of the quadratic and cubic action. As the non-attractor phase 
is associated with some paramteres at the background that evolves fast, one has to be careful 
in obtaining the actions through expanding these parameters. Note that the mechanism of 
breaking time-translation invariance is not different from usual attractor phase so we can use 
result of [3] prior to expansion which are presented in (]2]) and (|4]). 

First, we have to obtain the contributions from H(t + n) and H 2 (t + n) to first order 
of e, since in principle the time derivatives of the Hubble parameter generates powers of r) 
which can not be discarded because \r]\ ~ 1 during the non-attractor phase. These expansions 
are performed in (fTOjl and (fill) . Also other parameters in action may vary fast so we have 
to expand them and extract their leading behaviors. These terms are M A (t + 7r) which 
corresponds to speed of sound and c 3 (t + 7r) which is a free parameter of the theory not being 
fixed by the symmetries. The expansion of these two parameters are presented in (fT2j) and 

na. 

Equipped by expansions of all parameters in the action, the second order action is given 
as follows 

S 2 = j d 4 xa 3 M 2 {eH 2 ( ix 2 — 2W<9j7r — 27 t5N + 5N 2 —^<9j7n9j7r) + 2er]H 3 iT'k+ (62) 

SerjH^ix 2 + 6 eH 3 nSN — (1- -)eH 2 [k 2 + 5N 2 — 2tx8N] 

C 2 S 

+er] 2 H\ 2 - 2H8Nd i N i + - d.N^Nj] - 3 H 2 SN 2 } + 0(e 2 ) 
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Simplifying the above action yields 


s 2 = t + 


-f QHttSN — 2N l di7r —4’^» 7r ^» 7r ) 
c 2 a 2 

+ [ d 4 xa 3 M 2 H 2 { - 2SN8 i N i - 35N 2 } 


(63) 


As usual, A/V and IV* are non-dynamical in the sense that their equation involves no time 
derivatives. Solving their equations algebraically we obtain 


5N = eHn + 0(e 2 ) 

+ 0(e 2 ) (64) 

It is worth mentioning that these results are independent of the attractor or non-attractor 
assumption. Also the results are independent of the small or large variation of sound speed 
as s = c s /Hc s does not appear in the quadratic action. 

Having obtained 5N and N l one has to plug them back into action (1631) to find the final 
second order action in terms of n. As can be seen from Eq. (164)1 both 5N and N l are 
proportional to e. However, during the non-attractor phase e falls off exponentially so one 
expects that the contributions from 5N and N l in the action (ICuTl) to be sub-leading. Indeed, 
from (|S3D we see that the terms containing 5N and N l are higher orders in e compared to 
terms containing only n. Therefore, to leading order in power of e, one can safely neglect the 
contribution from 5N and N l in the action corresponding to set 5N = N l = 0. This represents 
our decoupling limit in which we keep only terms of 0(e) and discard terms of higher orders 
in e. Also in this limit, one can safely neglect the contribution from the gravitational action 
(]4)). In the decoupling action (163|) reduces to Eq. (TT4)l . 

The third order action is complicated in general but in the decoupling limit as described 
above the action simplifies significantly yielding 


& = 


d A xa 3 eM 2 H 2 


( T) — 2 s)H 


7T7T 2 — r]H — ll(di'Ky 


a 

- (i+^ 3 )(i - ^)tt 3 +(i - 

° C s c s a 


(65) 


Note that, unlike second order action, this action is sensitive to large time variation of c s as 
s appears directly in the action. In addition the cubic action depends on c 3 too. 


References 

[1] P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5076 [astro-ph.CO]. 

[2] P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5082 [astro-ph.CO]. 


19 








[3] C. Cheung, P. Creminclli, A. L. Fitzpatrick, J. Kaplan and L. Senatore, JHEP 0803, 014 
(2008) [arXiv:0709.0293 [hep-th]]. 

[4] L. Senatore and M. Zaldarriaga, JHEP 1204, 024 (2012) [arXiv: 1009.2093 [hep-th]]. 

[5] P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys. 571, A24 (2014) 
[arXiv:1303.5084 [astro-ph.CO]]. 

[6] M. Alishahiha, E. Silverstein and D. Tong, Phys. Rev. D 70, 123505 (2004) 

| hep-th /0404084j. 

[7] E. Komatsu et al. , “Non-Gaussianity as a Probe of the Physics of the Primordial Universe 
and the Astrophysics of the Low Redshift Universe,” arXiv:0902.4759 [astro-ph.CO]. 

[8] X. Chen, “Primordial Non-Gaussianities from Inflation Models,” Adv. Astron. 2010, 
638979 (2010) [arXiv:1002.1416 [astro-ph.CO]]. 

[9] E. Komatsu, “Hunting for Primordial Non-Gaussianity in the Cosmic Microwave Back¬ 
ground,” Class. Quant. Grav. 27, 124010 (2010) [arXiv:1003.6097 [astro-ph.CO]]. 

[10] J. M. Maldacena, JHEP 0305, 013 (2003) |astro-ph/0210603|. 

[11] P. Creminelli and M. Zaldarriaga, JCAP 0410, 006 (2004) |astro-ph/0407059|. 

[12] E. Pajer, F. Schmidt and M. Zaldarriaga, Phys. Rev. D 88, no. 8, 083502 (2013) 
[arXiv: 1305.0824 [astro-ph.CO]]. 

[13] M. H. Namjoo, H. Firouzjahi and M. Sasaki, Europhys. Lett. 101, 39001 (2013) 
[arXiv: 1210.3692 [astro-ph.CO]]. 

[14] X. Chen, H. Firouzjahi, M. H. Namjoo and M. Sasaki, Europhys. Lett. 102, 59001 (2013) 
[arXiv:1301.5699 [hep-th]]. 

[15] X. Chen, H. Firouzjahi, E. Komatsu, M. H. Namjoo and M. Sasaki, JCAP 1312, 039 
(2013) [arXiv: 1308.5341 [astro-ph.CO]]. 

[16] W. H. Kinney, Phys. Rev. D 72, 023515 (2005) |gr-qc/0503017 . 

[17] H. Motohashi, A. A. Starobinsky and J. Yokoyama, arXiv:1411.5021 [astro-ph.CO]. 

[18] C. Armendariz-Picon, T. Darnour and V. F. Mukhanov, Phys. Lett. B 458, 209 (1999) 

| hep-th /9904075]. 

[19] J. Garriga and V. F. Mukhanov, Phys. Lett. B 458, 219 (1999) |hep-th/9904176|. 

[20] Q. G. Huang and Y. Wang, JCAP 1306, 035 (2013) [arXiv: 1303.4526 [hep-th]]. 


20 


[21] C. Cheung, A. L. Fitzpatrick, J. Kaplan and L. Senatore, JCAP 0802, 021 (2008) 
[arXiv:0709.0295 [hep-th]]. 

[22] V. Assassi, D. Baumann and D. Green, JHEP 1302, 151 (2013) [arXiv: 1210.7792 [hep- 
th]]. 

[23] R. Flauger, D. Green and R. A. Porto, JCAP 1308, 032 (2013) [arXiv: 1303.1430 [hep- 
th]]. 

[24] L. Berezhiani, J. Khoury and J. Wang, JCAP 1406, 056 (2014) [arXiv:1401.7991 [hep- 
th]]. 

[25] L. Berezhiani and J. Khoury, JCAP 1402, 003 (2014) [arXiv: 1309.4461 [hep-th]]. 

[26] V. Sreenath, D. K. Hazra and L. Sriramkumar, arXiv: 1410.0252 [astro-ph.CO]. 

[27] V. Sreenath and L. Sriramkumar, JCAP 1410, no. 10, 021 (2014) [arXiv: 1406.1609 
[astro-ph.CO]]. 

[28] P. Creminelli, J. Norena and M. Simonovic, JCAP 1207, 052 (2012) [arXiv:1203.4595 
[hep-th]]. 

[29] P. Creminelli, A. Perko, L. Senatore, M. Simonovic and G. Trevisan, JCAP 1311, 015 
(2013) [arXiv: 1307.0503 [astro-ph.CO]]. 

[30] E. Dimastrogiovanni, M. Fasiello, D. Jeong and M. Kamionkowski, arXiv:1407.8204 
[astro-ph.CO]. 

[31] X. Chen and Y. Wang, JCAP 1004, 027 (2010) [arXiv:0911.3380 [hep-th]]. 

[32] X. Chen and Y. Wang, Phys. Rev. D 81, 063511 (2010) [arXiv:0909.0496 [astro-ph.CO]]. 

[33] R. Emami, JCAP 1404, 031 (2014) [arXiv:1311.0184 [hep-th]]. 

[34] S. Pi and M. Sasaki, JCAP 1210, 051 (2012) [arXiv: 1205.0161 [hep-th]]. 

[35] X. Chen, M. x. Huang, S. Kachru and G. Shin, JCAP 0701, 002 (2007) |hep-th/0605045|. 

[36] S. E. Shandera and S.-H. H. Tye, JCAP 0605, 007 (2006) |hep-th/06010991. 

[37] S. Weinberg, Phys. Rev. D 72, 043514 (2005) [hep-th/0506236]. 

[38] Y. Wang, arXiv:1303.1523 [hep-th]. 


21 


